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Abstract 

A simple position probability density formulation is presented for the motion of a particle in a 
spherically symmetric potential. The approach provides an alternative to Newtonian methods for 
presentation in an elementary course, and requires only elementary algebra and one tabulated in- 
tegral. The method is applied to compute the distributions for the Kepler-Coulomb and isotropic 
harmonic oscillator potentials. Formulas are also deduced for the average values for powers of 
the radial coordinate, and applied to describe perturbations to these systems. The classical re- 
sults are also compared with quantum mechanical calculations using the Einstein-Brillouin-Keller 
semiclassical quantization. 



I. INTRODUCTION 



A significant distinction exists between the conceptual framework presented in traditional 
introductory physics courses and that used in the advanced physics courses that follow 
themi. Introductory physics courses utilize historical Newtonian concepts involving forces 
and accelerations, but these concepts never enter in more advanced formulations. The 
introductory approach is often characterized as "classical" whereas that of the more advanced 
is described as "quantum mechanical." However, the primary difference between the two 
approaches arises not because of quantization, but instead from a nonessential heuristic 
tendency to describe macroscopic systems by instantaneous values for position, speed, and 
acceleration, and microscopic systems by time-averaged position probability densities. 

The reasons for this are clear, since a macroscopic trajectory is disturbed only slightly 
when successively interrogated with visible light, whereas a microscopic system may be 
destroyed by interrogation with a single short- wavelength photon. Thus the description of 
the microscopic system requires the superposition of many similarly interrogated systems. 
Unfortunately, this dichotomy produces a serious disconnect between physics as it is taught 
to non-major students in service courses and physics as it is practiced. Despite efforts to 
inject modern topics into a Newtonian presentation, this discontinuity further widens the 
gap between physics and society. 

In a recent essay, Wilczek^ has described the force concept as an insubstantial "culture" 
that provides a common language, but not an algorithm for constructing the mechanics 
of the world. Similarly, Taylor 3 has suggested an alternative approach that uses the least 
action principle in place of Newtonian forces. Both essays provide persuasive historical 
quotes from respected authorities who have urged that the force approach to the teaching 
of elementary physics be replaced. Unfortunately, the Newtonian model offers practical 
advantages, particularly in the testing and evaluation of student performance, and is thus 
very firmly entrenched. 

It is sometimes argued that initial use of the Newtonian approach is necessary, because 
a quantum mechanical formulation would be too demanding mathematically. However, the 
problems attacked in elementary textbooks tend to be simpler than those treated in quan- 
tum mechanical textbooks. If one examines problems of similar complexity, a Newtonian 
formulation is often much more complex mathematically than the corresponding quantum 
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mechanical solution. For example, elementary textbooks describe the two-dimensional Ke- 
pler orbit problem, but it is invariably restricted to the special case of a circular orbit (or, 
in the flat earth approximation, to a parabolic trajectory). When the classical problem 
is formulated in terms of position probability densities, three-dimensional elliptic orbits are 
automatically included. Moreover, deviations from a pure inverse square law can be included 
as perturbations 4 , all in a purely classical framework. It is also possible to add semiclassical 
quantization directly to the classical solution when desired. 

A formulation is presented here in which the periodic three-dimensional motion of a par- 
ticle in a central potential is treated in terms of classical position probability densities. The 
method is applied to the problems most frequently encountered in an introductory quantum 
mechanics course, namely the Kepler-Coulomb and isotropic harmonic oscillator potentials. 
While these two potentials lead to solutions that possess certain symmetries, they also have 
interesting differences. For example, the Kepler- Coulomb exemplifies an interaction that 
decreases with increasing separation, whereas the isotropic harmonic oscillator exemplifies 
an interaction that increases with increasing separation. 

In this presentation the position probability densities are evaluated, closed form expres- 
sions for the average values for powers of the radial coordinate are obtained, calculations are 
made for sample perturbations of the systems, and the connection to the EBK semiclassical 
quantization is prescribed. 

II. POSITION PROBABILITY DENSITIES FOR CENTRAL POTENTIALS 

Consider a particle of mass m moving in a central potential V(r) described by the standard 
spherical polar coordinates r, (p. For periodic motion with period T, the dwell time, or 
position probability density, is given by 



where p r is the radial component of the momentum of the particle, which can be described 
using conservation of energy as 




(1) 




(2) 
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Here L is the angular momentum. For a prescribed potential, the radial momentum can be 
obtained as 

Pr = V2m~\/E -V(r)- L 2 /2mr 2 . (3) 

With periodic orbital motion, the radial coordinate will undergo librations between turning 
points that are specified by the roots of the equation 

Er 2 - V(r)r 2 - L 2 /2m = . (4) 

For the potentials considered here there will be two roots to the equation, denoted as A±. 
Since the potential involves only r, the angular momentum will be constant over the orbit. 
In the case of the Kepler-Coulomb and isotropic harmonic oscillator potentials, the orbits 
are both ellipses, so Kepler's second law of equal areas swept out in equal times is valid for 
both. Thus 

l 2 d<p = 7rab = _L_ 

2 dt T 2m K } 

where a and b are the semimajor and semiminor axes of the ellipse, and nab is its cross 
sectional area. This equation permits the specification of the period, which provides the 
normalization of the distribution. If N denotes the number of librations in a period (N=2 
for the Kepler-Coulomb, N=4 for the harmonic oscillator), then the average values of powers 
of r are given by 

AT r A + 

(r k ) = — / drP(r)r k (6) 



A. Kepler- Coulomb Potential 

The potential 

V{r) = -k/r (7) 

gives rise to a negative (binding) energy, which we denote as E B = —E so as to explicitly 
display the sign within square roots. The momentum thus becomes 

p r = V2m^E + k/r - L 2 /2mr 2 (8) 

with turning points given by the roots of 

- E B r 2 + kr- L 2 /2m = (9) 
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given by 

h ( b \ 2 r.2 

(10) 



A —-<- ( k V 12 
± 2E B \ \2E B ) 2mE L 



In this case the coordinate system is centered on one of the foci of the ellipse, for which the 
semimajor and semiminor axes are given by 

a = k/2E B 



b = L /^2mE B . (11) 

An example of such an orbit with a = 1 unit and b = a/2 is shown in Fig. la. 
The period can be computed from the definition of b using Eq. El in the form 



m 2mnab 2m , , 

T = — = ™Vi£- (12) 

Inserting these relationships into Eq. El (with iV=2 since here the periapsis and apoapsis are 
separated by 180°) 

I Tflv 

P(r)dr = — , . (13) 

™ y/{A+-r){r-A-) 

The position probability density corresponding to the orbit in Fig. la is shown in Fig. 2a. 



B. Isotropic harmonic oscillator 

The potential 

V(r) = kr 2 /2 (14) 

yields the momentum 



p r = V2myJE - kr 2 /2 - L 2 /2mr 2 . (15) 
with turning points specified by the roots of the equation 

Er 2 — kr A /2 — L 2 /2m = , (16) 

given by 

^ = W(i) (17) 

This orbit is also elliptical, and is comparable to that of the Kepler- Coulomb system, except 
for the fact that the coordinate system is at the center of the ellipse rather than at one of 
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the foci. Here the turning points are at the semimajor and semiminor axes 



a = A + 

b = A. , (18) 

and the period corresponds to four of these turning points. An example of such an orbit, 
also with a = 1 unit and b = a/2, is shown in Fig. lb. 
The area of this ellipse is 



nab = Ja 2 + AL = . (19) 

ymk 

Using Eq. E3 this gives a value for the period 

T= 2 -^ = 2n^. (20) 

Inserting these relationships into Eq. El (with N=4 since here the closest approach and 
furthest recession are along the semiaxes, and thus separated by 90°) 

P(r)dr = - ; ^ . (21) 

* J(Al-r 2 )(r 2 -A 2 _) 

The position probability density corresponding to the orbit in Fig. lb is shown in Fig. 2b. 



III. EXPECTATION VALUES 

Average values of quantities weighted by these distributions can be obtained by directly 
integrating these expressions. However, they can also transformed into the form of the 
standard integral 5 

1 r 2 <P ,1 

— / Mi+,cos^r = (i-e 2 r/ 2 p n (-==) (22) 

27r Jo Vl — e 2 

where P n (x) is the Legendre polynomial (in an unusual application where the argument 
x > 1). Negative powers can be handled using the relationship 

P_ n (x) = P n _ x {x) . (23) 

In addition to the radial integral formulation of Eq. H3 the expectation value can alterna- 
tively be written as 

( r *) = I ( T dt r k = 1 r k . (24) 

w T Jo T Jo d<p/dt v ; 
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Conservation of angular momentum relates r and ip through Eq. which can be rewritten 

T d<p/dt = 2nab/r 2 . (25) 

Inserting this into Eq. |2H 

rt=ir^^ a - (26) 



27ra6 jo 

It remains only to choose the equation of the orbit, and to use Eq. 1221 to evaluate this 
expectation value. 

A. Kepler- Coulomb problem 

Here the coordinate system is centered on one of the foci of the ellipse, which has the 
equation 

~ = ^(l + £cos^) (27) 



where e = J 1 — b 2 /a 2 is the eccentricity of the ellipse. Inserting this relationship for r into 
Eq.HU 



1 f a\- k - 2 1 r 27T 

( 

which, using Eq. [221 becomes 



^ - bt (!) p - (?) • < 29 > 



A few examples are: 



(r) 


= a [3 - (6/ 




= 1/a 




= 1/ab 


(r" 3 ) 


= 1/b 3 


(r- 4 ) 


= (r)/b 5 . 



(30) 



B. Isotropic harmonic oscillator problem 

In this case the center of the coordinate is at the center of the ellipse, and has the 
equation^ 

1 1/1 1\ /l 1\ 

V 2 = 2^ 2 + b 2 )-\^-b 2 ) C0 ^ (31) 
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which can be rewritten 



r 



£ / g 2 + ^ 

2 ~ I 4a 2 6 2 



Defining here 



i 2 -b 2 



a 2 + b 2 

the quantity occurring in Eq. |221 simplifies to 



\ 



1 



'a 2 -b 2 ' 
a 2 + b 2 



2ab 
a 2 + b 2 



The expectation value is given by 



(n = -i 



1 (a 2 + b' 



k + 2 

2 i h2\ — 



1 r 2 * 
— / dip (1 + e cos 2<f) 
2n Jo 

a 2 + b 2 \ 



ab \ 2a 2 b 2 



which integrates to 



(r k ) = (ab) k/2 P_M2 ( 



k + 2 



2ab 



(32) 



(33) 



(34) 



(35) 



(36) 



This result is valid for both odd and even powers. For odd powers, the Legendre function 
can be evaluated numerically as a hypergeometric series, as shown in the Appendix. 
A few examples are: 



(r 2 ) 
(r 4 ) 
(r- 2 ) 
(r- 4 ) 



a 2 + b 2 )/2 

3 (a 2 + b 2 ) - Aa 2 b 2 } /i 
l/ab 

(a 2 + b 2 )/2a 3 b 3 . 



(37) 



IV. PERTURBATION CALCULATIONS 

One of the strengths of this method is the ease with which perturbations to the energy 
of the system can be computed. The total energy can be deduced from the potential using 
the virial theorem 

E = (V(r)) + \ (r^pj (38) 
so a perturbation of the form AV(r) can be computed as 



E' = E+(AV(r)) 



(39) 
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A. Example 1: Kepler- Coulomb with a 1/r 3 perturbation 



This can occur, for example, in an atom with a spin-orbit magnetic interaction, or in a 
gravitational system with a Schwarzschild general relativistic correction^. 
The energy of the system is 

E = (-kr- 1 ) + ^(kr- 1 ) (40) 
If the perturbation is AV(r) = A/r 3 , the perturbed energy is 

E> = - k ~(r-') + A(r- 3 > 

= "25 + 65 (41) 
which results in a precession of the ellipse. 

B. Example 2: Anharmonic oscillator with an r 4 perturbation 

The energy of the system is 

E = (hr 2 ) + (hr 2 ) (42) 
If the perturbation is AV(r) = Ar 4 , the perturbed energy is 
E' = k(r 2 ) + X(r 4 ) 



2 V ^ J ^ 4 
which also results in a precession of the ellipse. 



3(a 2 + 6 2 ) 2 -2a 2 6 2 | (43) 



V. THE SEMICLASSICAL EBK QUANTIZATION 

The semiclassical Einstein-Brillouin-Keller quantization is given by 

K + = j dq { pi (44) 

where \x the Maslov index, which is the number of turning points. This formalism was 
applied for spherical symmetric potentials in an earlier paper—. The angular phase integrals 
yield a value for the angular momentum 

L=(£ + l/2)h. (45) 
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The square of this result 

L 2 = [£(£ + 1) + 1/4] h 2 (46) 

agrees with the quantum mechanical result in the correspondence limit. 

Our earlier calculations 8 for the radial phase integral permit the specification of the 
semimajor and semiminor axes of the ellipses. 



A. Kepler- Coulomb 

For the Coulombic atomic problem (k = Ze 2 /47reo for a hydrogenlike atom), 



(47) 



The radial quantum number n r is displayed here so that the two potentials can be com- 
pared under conditions whereby n T and £ have the same range of values 0, 1, 2, ... The 
expression is usually written in terms of the principal quantum number n = n r + £ + 1. 

With this quantization the perturbed energy of Eq. PHI becomes 



E' 



k 2 m 



1 



2A 



1 



k n 3 (l 



1)2 

2/ 



(48) 



which agrees with the quantum mechanical result with the correspondence (£ + 
£{£ + \){£ + l) . 



1^3 



B. Isotropic harmonic oscillator 



In this case the quantization yields value for the semiaxes (denoting u = yk/m 

a 2 + b 2 h ( „ 3 
2n r + £ + - 



2 moo \ 2, 

ab = —(£+-) . (49) 

Here again the radial quantum number n r is displayed for comparison with n r and £ 
having the same range of values 0, 1, 2, ... The expression is usually written in terms of the 
quantum number n = 2n r + £ . 
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With this quantization the perturbed energy of Eq. H3] becomes 



E' = %u [n + - 



3\ xn 2 



3\ 2 /„ V 2 



3(n + -j 



(50) 



2) 2m 2 uj 2 

which agrees with the quantum mechanical result^ with the correspondence [I + |) 2 — > 
(£-l)(£ + l). 

VI. CONCLUSION 

This formulation in terms of the classical position probability density provides a mathe- 
matically simple exposition of the difference in frameworks between classical and quantum 
mechanical physics. Although this one exercise does not provide a comprehensive alternative 
to the standard presentation, it can clearly illustrate at the introductory level the limitations 
of the Newtonian approach. 

APPENDIX A: LEGENDRE FUNCTIONS OF HALF-ODD-INTEGER ORDER 

Legendre functions of half-odd-integer order can be evaluated using the hypergeometric 
series 

P-*-i(z) = PAz) = (i±£)" F(-u,-v; 1; ^j) . (Al) 



Thus 



a + b\ k „/ k k ( a — 6 X 2 ^ 



rt = m * ur* •■ (A2) 



For the case shown in the figures, b = a/2, this gives for the first moment, 

3a / 1 1 1 \ . t . 

(r) = — r F [—-,—- ; 1; - ) = 0.77098 a . (A3) 



4 V 2' 2 ' ' 9 

In the limit b —* we can use the fact that F(a, a, ; c; 1) = r(c)T(c — 2a)/T(c — a) 2 to write 
the moments (k > 0) for a linear oscillator in one dimension: 

k\ cd^ 

{rk)l ° = axTTW ' (A4) 

These results check against the elementary results, for example 

(r) 1D = 2a/n 
(r 2 ) 1D = a 2 12 

(r 4 ) 1D = 3a 4 /8 . (A5) 
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FIGURE CAPTIONS 

Figure 1. Comparison of the elliptic orbits with a = 1 unit and b = 0.5 for the two 
examples. 

Figure 2. Classical position probability distributions for the two elliptic orbits shown in 
Fig. 1. 
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